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^— ^ ■ Abstract: We compute the mean free path and shear viscosity in the color-flavor locked 

■^ ' (CFL) phase of dense quark matter at low temperature T, when the contributions of 

^""^ mesons, quarks and gluons to the transport coefficients are Boltzmann suppressed. CFL 



quark matter displays superfluid properties, and transport phenomena in such cold regime 
are dominated by phonon-phonon scattering. We study superfluid phonons within thermal 



qj I field theory and compute the mean free path associated to their most relevant collision pro- 

'^ ' cesses. Small-angle processes turn out to be more efficient in slowing transport phenomena 

in the CFL matter, while the mean free path relevant for the shear viscosity is less sensi- 

^ ■ five to coUinear scattering due to the presence of zero modes in the Boltzmann equation. 

j^ ■ In analogy with superfluid He^, we find the same T power law for the superfiuid phonon 

damping rate and mean free path. Our results are relevant for the study of rotational 

properties of compact stars, and correct wrong estimates existing in the literature. 
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1. Introduction 

It has been known for long time that cold dense quark matter should exhibit the phe- 
nomenon of color superconductivity, at least at asymptotic large baryonic densities [1]. 
Only recently this interesting scenario has been considered thoroughly, both with the hope 
of deepening our understanding of QCD, and for its applications in astrophysics (see [2, 3] 
for reviews and references). These recent studies realized that the fermionic gap A might 
be as large as 100 MeV, for quark chemical potentials /i of the order of 500 MeV. Thus the 
phenomenon should lead to clear distinct features at the macroscopic level. Much efforts 
are now being devoted to finding signatures of color superconducting quark matter in any 
of its different possible phases. There are different superconducting phases, according to 
the number of quark flavors that participate in the diquark condensation. In this article 
we will only be concerned about the color-flavor locked (CFL) phase [4], which occurs in 
the presence of three light quark flavors. 

An important set of signatures of color superconductivity might be found in phenom- 
ena associated to stellar vibration and/or rotation. The existence of r(otational)-mode 
instabilities in all relativistic rotating stars [5] seems to be incompatible with the existence 



of millisecond pulsars, unless the instabilities are suppressed by sufficiently large viscosities. 
Thus, the values of the viscosities can be used to rule out models for millisecond pulsars. 
While the common belief is that pulsars are neutron stars, it cannot be discarded that 
they are quark or hybrid stars. Madsen did point out that millisecond pulsars could be 
made of strange quark matter [6], characterized by large shear and bulk viscosities, but 
he ruled out [7] quark matter in the CFL phase. This last conclusion, however, was based 
on a wrong estimate, as it was assumed that the main contribution to the viscosities was 
due to gapped quarks, being exponentially suppressed as exp[— 2A/T], for values of the 
temperature T of the order of 0.1 MeV and below. 

Transport properties in the CFL phase of QCD are not dominated by the quarks. In 
this phase the diquark condensate breaks spontaneously the baryon symmetry [/^(l), and 
CFL quark matter is a superfluid. Chiral symmetry is also spontaneously broken. Asso- 
ciated to those breakings, there are Goldstone bosons, which are light degrees of freedom. 
In addition, there is an unbroken U{1) subgroup whose gauge boson, a combination of the 
photon and one gluon, is massless at zero temperature and can be viewed as the in-medium 
photon. A CFL quark star is electrically neutral [8], both for the real and in-medium elec- 
tromagnetism, but at finite temperature one may also expect to find electrons. All the 
above mentioned particles are light and their contribution to the transport coefficients in 
this phase is bigger than that of the gapped quarks. Let us mention that we are considering 
a CFL quark star after the deleptonization era, so that all the neutrinos have escaped from 
the star. At higher temperatures, the contribution of neutrinos to the transport coefficients 
could be estimated from their mean free path, as computed in refs. [9, 10], but we will not 
consider such a temperature regime here. 

Chiral symmetry is not an exact symmetry of QCD. Therefore, the associated (pseudo) 
Goldstone bosons (the pions, kaons, and etas) are massive [11, 12]. At asymptotic large 
densities, meson condensation may occur [13, 14]. At more moderate densities, instanton 
effects become relevant, and modify the meson mass pattern [15, 16]. There are still some 
uncertainties in the instanton contribution to the meson masses, however some primary 
computations indicate that all meson masses m.,^ are in the range of tens of MeV, or even 
larger [16]. One can safely state then that at very low temperatures, T <C 1 MeV, as for 
the estimates of ref . [7] , the contribution of the mesons to the transport coefficients is also 
very suppressed by the Boltzmann factor exp[— m^/T]. 

In the regime T <C A, m.,^, light particles dominate transport. This was already pointed 
out in refs. [17, 18], where estimates for the thermal and electrical conductivities were 
given. The contribution of the in-medium electromagnetism to the shear viscosity can be 
extracted from that of a QED plasma, with minor modifications, as we carefully explain 
in Sect. 6, but it is small. We report here a new computation of the contribution to the 
shear viscosity of the only truly massless Goldstone boson of the CFL phase, the superfluid 
phonon. Its self-interactions have been derived in an elegant way by Son [19]. 

Shear viscosity describes the relaxation of the momentum components perpendicular 
to the direction of transport, and it is usually dominated by large-angle collisions. In 
this manuscript we compute the mean free path associated to both small and large angle 
collisions. The differential cross section of binary collisions mediated by phonon exchange 



is divergent for small-angle collisions. This is the typical Coulomb- Rutherford collinear 
divergence induced by massless exchange. In an ordinary scalar theory such a divergence 
does not appear, as a thermal mass is generated even if the boson is massless in vacuum. 
But the phonon remains massless at finite temperature, as thermal effects do not represent 
an explicit breaking of baryon symmetry. The divergence is regulated by the finite width of 
the phonon, or more precisely, by Landau damping, a process only occurring in a thermal 
bath. After regularization, we find that small-angle processes have a shorter mean free 
path than large-angle ones. This suggests that they might be more relevant for transport, 
as a large-angle collision can always be achieved by the addition of many small-angle ones. 
To compute the shear viscosity we solve the Boltzmann equation, linearizing it for small 
departures from equilibrium. Then one finds that a zero mode occurs in the collision 
operator, which suppresses severely the contribution of exact collinear processes. Then one 
needs to extend the integration over phase space beyond the strictly collinear limit. The 
full analysis is rather technical, and requires numerical evaluation of all processes. We also 
find that the phonon of the CFL and He^ superfluids share many properties, and have 
the same power laws for their damping and mean free path, suggesting a sort of universal 
behavior. 

This article is structured as follows. We review the superfluid hydrodynamic equations 
and the phonon effective Lagrangian in Sec. 2. In Sec. 3 we compute the phonon self-energy 
at one-loop. From the imaginary part of the on-shell energy we evaluate the phonon 
damping rate. In Sec. 4 we present the mean free path of splitting processes and binary 
collisions. The computation of the shear viscosity is done in Sec. 5, and in Sec. 6 we explain 
why the contribution of the in-medium electromagnetism is negligible at low temperatures. 
We comment in Sec. 7 on possible corrections to our results, and present conclusions in 
Sec. 8. We leave for the Appendices A and B the technical and numerical details of several 
computations. 

2. The CFL superfluid and the phonon low-energy eff'ective theory. 

The low energy effective theory for the only truly Goldstone boson of the CFL phase can 
be constructed from the equation of state (EOS) of normal quark matter [19]. For three 
massless quark flavors the EOS is 

Pif^) = ^/ , (2.1) 

where ^ is the quark chemical potential. From Eq. (2.1) Son obtained the effective La- 
grangian for the superfluid phonon 

CeS = ^[{do^-^f-{d,^f]' . (2.2) 

There is an interesting interpretation of the equations of motion associated to f, as 
they can be re-written as the hydrodynamical conservation law of a current representing 
baryon number flow, 

d^^imu'') = , (2.3) 



where no = -d7[\i^=fj.o is interpreted as the baryon density [19]. The superfluid velocity u'^ is 
proportional to the gradient of the condensate phase, as in Landau's model of superfluidity 

[25], 

u. = -^, (2.4) 

where D,yip = dyip — (^,0), and /io = {D^ipD'^ip)^''^. The energy momentum tensor can 
also be written in terms of the velocity defined in Eq. (2.4) and Noether's energy-density 

T^P = (e + P)u''uP - g^fP . (2.5) 

It is conserved and traceless 

d^T^'P = , T; = . (2.6) 

Eqs. (2.3) and (2.5) are the hydrodynamic equations for the relativistic superfluid 
[19]. They need modifications at finite temperature T, as phonons are thermally excited 
and conform a different fluid in the system [26, 27]. The well-known two-fluid model is 
necessary, featuring a superfluid that flows perfectly and shows no dissipation, and a normal 
fluid where dissipative processes are possible. At low T, these will be controlled by the 
thermal properties of the phonons composing the normal fluid and their scattering rates, 
that we compute in the following sections. To proceed with those calculations we first 
rescale the phonon field 

<P=^-^^ (2.7) 

vr 

to normalize the kinetic term in accordance with the LSZ formula. Then the Lagrangian 

for the rescaled field reads 

C,s = \{d,(t>? - y (5.0)=^ - ^doct>{d,<pdP<t>) + Y^(^M'^^'''^)' ' (2-8) 

where v = l/v3 is the phonon velocity. We have neglected above a total time derivative 
term, which is irrelevant for all computations to follow. The Lagrangian (2.8) is not 
renormalizable, and only valid up to scales of the order 2A, the energy necessary to break 
a diquark condensate. Eq. (2.8) only contains the lowest derivative couplings, valid for low 
energy scales. Zarembo [20] has shown that to all orders in — , the interaction terms are 
suppressed by powers of /i. Thus, for the temperature regime we will consider here T ^ A, 
those would only give very tiny corrections. 

Also of interest to us is Zarembo's observation [20] that the phonon dispersion relation 
at low momenta, beyond Son's theory, becomes 



^s{k) = -7=k 



11 k' .^fk''\ 



^ 540A2+^VaV 



(2.9) 



As k is increased, the phonons move slower, with the tendency of suppressing collinear 
splitting (a phonon cannot decay into two phonons of larger joint energy). We will come 
back to this point in the following sections. 



n = 




Figure 1: One- loop contributions to the phonon self-energy. 



3. One-loop phonon self-energy 

In this Section we compute the phonon self-energy, and evaluate its damping rate at one- 
loop. We use the Imaginary Time Formalism (ITF) to perform the calculation. Feynman 
rules are easily deduced from the Lagrangian (2.8). 

The phonon propagator in the ITF with momentum K = (A;o,k) = (2a;„,k), reads 

S{K) ^ -^ , (3.1) 

where w„ = lirnT, with n G ^, is a bosonic Matsubara frequency, and E^ = vk, where 
A; = |k|. 

There are two different diagrams that contribute to the one-loop self-energy, see Fig. 
1. For external momentum P = {po,p) = {iuj,p), these are 

and 

respectively, where 

F{P,K) = {po{2K-P-K^) + ko{P^-2K-P)f . (3.4) 

Above we have used the notational convention of the ITF 



/ 



^ 2^ / 7^rv3 ' (3-5) 



(2ttY ~ ^ J (2^)3 ' 



where the sum is over the Matsubara frequencies. We use dimensional regularization to deal 
with the ultraviolet (UV) divergences of the T = part of the diagrams, so that d = 3 — 2e. 
Since the thermal part of the diagrams is UV finite, we will analytically continue back to 
d = 3 for their evaluation. 

The superfluid phonon is a Goldstone boson, and since thermal effects do not represent 
an explicit breaking of the C^b(1) symmetry, its self-energy should vanish at P = 0. It is 
actually easy to check that this holds at one-loop level 

n(") (p = 0) = n^'') (p = 0) = , (3.6) 



so no thermal mass is generated. This property of the self-energy should hold to all orders 
in perturbation theory. 

After performing the sum over Matsubara frequencies we find that the first diagram 
only corrects the phonon velocity by a term proportional to (T/^)^ <^ 1. 

We proceed with the evaluation of the sum of Matsubara frequencies for the second 
diagram. We find 






Sl,S2=± 



S1S2 l + f{siEi)+f{s2E', 



1, 



ko=siEi \4E1E2 iuo — SiEi — S2E2 

(3.7) 
where Ei = vk and E2 = v\p — k|, and 

m ^ -^^ . (3.8) 

After analytical continuation to Minkowski space with retarded boundary conditions, 
iuj ^ Pq + ie, one notes that this diagram has both real and imaginary parts, the last one 
being related to the damping of the phonon. 

The real part of the diagram gives contributions that behave as follows for low mo- 
menta. For pq = p, the thermal corrections are proportional to {T/fj,)'^p'^, while the T = 
corrections go as p^/n In (M^/p^), where M is a renormalization scale. As for the previous 
diagram, these corrections are very much suppressed as compared to the tree level physics 
for p <^ T <^ fi, and we will neglect them. 

Although we use the full imaginary part of the self-energy in our computer work, 
numerically evaluating Eq. (3.7), we can display an exact analytical limit for illustration. 
For low external momentum, po,p <^ T, we keep only quadratic terms in the external 
momentum in F, so that 

lrnnip„p)^^J^^^l^{il + f, + f2){sipo-E,-E2) (3.9) 

- 6{po + Ei + E2)) - (/i - /2) (<5(po -Ei + E2) - d{po + Ei- E2)) } , 

where /j = f{Ei). For p^^p ^ it is easy to check that the last two delta functions provide 
the leading order behavior to the integral, so we neglect the first two deltas. For the last 
two delta functions, we perform the following approximations 

£'1 - £'2 ~ f p • k , /i - /2 ~ !> P • k--^ , (3.10) 

aEi 

where k = k/A;. Thus, one finds for d = 2> 



and after evaluating the integrals, one gets 

imn(po,p)-^:l^4e(.V-pg) , (3.12) 



where is the step function. This imaginary part of the self-energy corresponds to what 
is known as Landau damping. Particles disappear or are created through scattering in the 
thermal bath, and not via the processes which occur at zero temperature. 

The damping rate is defined in terms of the on-shell imaginary part of the self-energy 



v^ 



l{Ep) = -^^lmU{po = Ep,p) . (3.13) 

Since we neglect the real corrections to the self-energy, we can take pQ = Ep = p/\/3. In 
the limit p — > 0, we find 

7=-— 7=— P + 0(p3). (3.14) 

405V3 /U* 

Let us insist that the T = contribution to the damping rate is subleading in p. The optical 

theorem indicates that this should be proportional to ~ p^ /fi'^, as the on-shell imaginary 

part of the self-energy is related to the square of the tree-level scattering amplitude. Thus, 

at very low momenta, the thermal damping dominates over the T = one (but again, we 

include all effects in our numerical analysis). 

Remarkably, we find that the phonons of the CFL and He^ superfluids share many 

properties. They both travel at the speed of sound of the system. Their damping rates 

for low momenta and at low T follow the same law, going as ~ pT'^ [22]. Their damping 

rate at T = also behaves as ~ p^ in both cases [23]. These universal properties are more 

easily spelled out if one writes an effective field theory for the Goldstone mode of the non- 

relativistic superfluid [24], where one finds the same sort of cubic and quartic derivative 

couplings as for the CFL superfluid. This paralellism will be further studied in a separate 

publication. 

4. Phonon mean free path 

In this Section we present explicit computations of the mean free path associated to 1 ^^ 2 
and 2 •i-^ 2 phonon collisions. This is a necessary pre-analysis for the study of transport 
coefficients in the system. In a dilute gas, viscous and other transport coefficients are 
proportional to the mean free path, and hence inversely proportional to the damping rate or 
scattering cross section between quasiparticles. Here diluteness implies that the mean flight 
time between collisions is much larger than the duration of the collision, and in addition 
Boltzmann's molecular chaos hypothesis (no correlation between successive collisions) must 
hold. 

4.1 Collinear splitting processes 

We evaluate the mean free path associated to splitting collinear processes of 1 ^^ 2 particles. 
These processes are kinematically forbidden for massive particles, but they are not for 
massless ones. Because the phonon dispersion relation is linear with Son's lagrangian, 
these processes are perfectly collinear, that is, the momenta of the three particles are 
aligned. Thus, they are not efficient for shear viscosity. They could be efficient if the 
phonon dispersion law was changed so that the process were not perfectly collinear, but 
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Figure 2: Tree level Feynman diagrams contributing to 2 ^^ 2 phonon scattering processes. 



as discussed in Sec. 2, corrections to Son's lagrangian stemming from higher derivative 
interaction terms in the effective expansion (suppressed in the T/// counting), disallow 
non-collinear splitting as shown by Zarembo [20]. 

The inverse of the mean collision time for these processes is obtained averaging the 
one-loop damping rate computed in Sec. 3. Thus 



1 

1^2 — — 

n 



where 



n 



d'^p 



(fp 
(2^ 



fiEp 



l{Ep)fiEp 



T3 



■C(3) , 



(4.1) 



(4.2) 



(2^)3-'^ P' 4^V 

and C{x) is the Riemann <^ function with ^(3) = 1.202.... 

The mean free path is given by A = ^. From simple dimensional analysis, it can be 
inferred that the mean free path for these processes increases at low T as fi^/T^. We 
have numerically computed T and the mean free path with the exact one-loop self-energy, 
checking this power law for low T (see Fig. 3). This may fail at higher T as perturbative 
corrections to Re 11 turn relevant, but that temperature regime is not considered in this 
article. 

4.2 Two-body elastic scattering processes 

Here we compute the mean free path of binary collisions, 2 ^^ 2 processes as depicted in 
Fig. 2. To evaluate their mean free path one needs the imaginary part of two-loop self- 
energy diagrams. We employ instead the more direct formulation of kinetic theory, which 
provides the same answer. Thus, the inverse of the mean collision time is given by 



1 2— »2 — 

2n 



where 



/ /p/k (1 + /p') (1 + /k') (2vr)45(4) iP + K-P'- K') |T|2 , (4.3) 

Jp.k.p'.k' 



Sp 



p , 2Ep{2TTf 



(4.4) 



and |Tp is the scattering amplitude squared. We have also used as shorthands /p = f{Ep), 
etc. 



Binary collisions with large or small-angle deflection behave very differently. Let us 
start with the first ones. Dimensional analysis simply suggests that these proceed with a 
mean free path of order fi^/T^ [17]. For the typical values of // ~ 500 MeV and T ~ 1 
MeV, the mean free path is then larger than 10^ km, and exceeds the typical radius of 
a compact star, i? ~ 10 km. This means that these processes are totally irrelevant for 
transport phenomena inside the star. 

On the other hand, let us consider the t-channel scattering in a binary collision, see 
Fig. 2 (the n-channel has the same behavior). It is easy to check that the associated scat- 
tering matrix diverges for small-angle collisions. In a massless scalar theory this divergence 
would be cured through the generation of a thermal mass in the medium, which corrects 
the low momenta behavior of the propagators. This cannot happen for the phonons de- 
fined in Sec. 2, as these particles remain massless at finite T. However, the phonon width 
regulates the divergence, as we detail below. 

Let Q = (a;, q) be the momentum transfer, Q = P — P' . Then the scattering matrix 
in the t-channel reads 

1^,2 Mtt^V F{P,Q)F{K,-Q) 

where F was defined in Eq. (3.4). We have only included the imaginary part of the phonon 
self-energy, as the real part only gives subleading corrections as already discussed. Were 
Imll = 0, the scattering matrix would diverge when uj"^ -^ f^g^, as the exchanged phonon 
would be on-shell. This corresponds to a collision where the momenta of the scattered 
particles is only slightly deflected by the collision. This divergence is regulated by including 
the damping of the exchanged phonon, which effectively amounts to a resummation of a 
certain class of diagrams. Because the particles in the thermal bath have typical energies 
of order T, the momentum transfer in the collision will behave typically as g <C T. So for a 
leading-order estimate of the mean free path, it would suffice to include the value given in 
Eq. (3.12), although for numerical computations we can calculate with the full self-energy, 
and check our approximations. 

After regularization, the 2^2 damping rate becomes finite, and dimensional analysis 
of the mean free path concurs with the full calculation on a fi^ /T^ power law (see Ap- 
pendix A). The numerical estimate (see Fig. 3) teaches us that the mean free path of this 
process is only slightly smaller than that of the splitting collinear processes. 

We have excluded from the analysis the s-channel, as it gives a subleading contribution 
to the mean free path. As opposed to the t and u channels, where the momenta transferred 
can be arbitrarily small, here it is of order T. Thus the regulated divergence in the s-channel 
has much less phase space in the integration region of Eq. (4.3), and can be safely put aside. 

Again, we have found the same T dependence for the phonon mean free path as 
in supcrfluid Helium below 0.6 K [21]. Initially, Landau and Khalatnikov assumed that 
the viscosity mean free path was governed by four phonon processes, scaling as T~^. 
Experimental values of the viscosity showed that this was wrong, and that the mean free 
path was dominated by small-angle collisions, behaving instead as ~ T^^, as in the CFL 
superfiuid. 
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Figure 3: Mean free path below T = 1 MeV, scaling as ^'^/T^. At very low T it exceeds the 
typical radius of a compact star. A hydrodynamical description in this regime is not meaningful 
anymore. 



5. Shear viscosity from phonon collisions 

In this Section we compute the shear viscosity in the phonon fluid employing kinetic theory. 
The shear viscosity rj enters as a dissipative term in the energy-momentum tensor as follows 



v{diyj-djVi--V -Y 



(5.1) 



where V is the three dimensional velocity of the phonon fluid in the frame where the 
superfluid component is at rest. The phonon four velocity is given by U^ = (7, —7V), with 

The phonon distribution function /, under the hypothesis of molecular chaos, obeys 
the transport equation 

(5.2) 



f-^[^^] 



In the frame where the superfluid is locally at rest, the Liouville operator appearing in 
Eq. (5.2) (and thus the equilibrium solutions) depends on the superfluid velocity Eq. (2.4). 
Because we are interested in the shear viscosity of the phonon fluid, and we will linearize 
the transport equation in gradients of V, we will neglect the superfluid velocity in all the 
developments to follow. We note, however, that this approximation would not be valid for 
the computation of other transport coefficients. 

The collision term refers to the binary collisions and splitting/joining processes de- 
scribed in Sec. 4. In the shear viscosity calculation, we neglect the last one for the reasons 
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explained in Sec. 4.1. The binary collision term is given by: 

C[U] = \I {27r)'6^'\P + K-P'-K')\T\'-LD, (5.3) 

with 

D = /p'/k'(l + /p)(l + /k) - /p/k(l + /p')(l + /k') . (5.4) 

As usual, the collision term vanishes for the equilibrium distribution /*"^', which is 
given by 

where P^ = (-Ep,?). 

For the computation of transport coefficients we have to consider small departures 
from equilibrium so that the distribution function can be written as 

f = n + 6f. (5.6) 

Accordingly, D can be written as D = D^"^' + 6D. 

To solve the transport equation we employ the Enskog expansion. To first order 
we consider the perturbation only on the collision term and linearize it, whereas in the 
advective term we take the equilibrium distribution function. Since the equilibrium part 
f^^' annihilates the collision term, we substitute D in Eq. (5.3) by 5D, that can be written 
as 

SD = /^^•/^<i-/;?-/^?-e(^^+^'=)/^A[-^(l - e-^/^)] , (5.7) 

where the symbol Ag is a shorthand for: 

Ag^g{p')+g{k')-gip)-g{k) . (5.8) 

From the kinetic point of view, the dissipative part of the spatial momentum-stress 
tensor can be written as: 

^fW/'P^ 5/ • (5.9) 

To simplify the computations, we assume a velocity profile of the form V = (V^(y), 0, 0). 
We then assume that the driving shear departures from the equilibrium distribution func- 
tion can be parameterized as 

Sf = 7^9{p)PxPy-^ ■ (5.10) 

Solving the linearized transport equation for g{p) would yield the shear viscosity, which 
is expressed as 

^ J^ f _JP5^^ feci. . (5.11) 

This expression can be viewed as a scalar product 

V = {xm (5.12) 
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with 

X = PiPj (5.13) 

and $ a solution to the Boltzmann equation that can be written as 

Ix) = C|<l>) (5.14) 

with C the hnearized colHsion operator. Equivalent to Eq. (5.12) is then 

ry=($|C|$) (5.15) 

that squared, and substituting (5.12), yields 

In this expression, <I> is the exact solution of the Boltzmann equation. But it is also 
amenable to a variational treatment considering a set of test functions $, where the exact 
solution $ would make (5.16) reach an extremum (as can be seen by substituting in it 
$ = $+e and observing that the linear terms in e cancel if $ solves the transport equation). 
In a rotationally invariant way, we can take the scalar product necessary in these 
expressions as 



E T W 



and we take as trial function 



^ = 9ij{p) = 9{p)Pij = 9(.P){PiPj - 3<^ij/) (5-18) 

dip) - 



1 _ e-Ep/T 



with 1/ a variational parameter. This family is general enough to allow treatment of a 
polynomial family such as used in [28], but also includes rational functions. For this 
family, Eq. (5.16) reduces to 



with 



and 



VM = ^ , (5.19) 



L[u] = -^ [ {27r)^ 6\P + K - P' - K')^-^ (5.21) 



40T yp,k,p',k' 2 

As the integrand of L[iy] is positive definite, so is also Eq. (5.19). The computation of /[z^] 
is easy, but to compute L[iy] requires a multidimensional Montecarlo integration. This is 
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quite a technical calculation and is detailed in the Appendices, but we will point out some 
subtleties here. 

It is worth stressing that for a quick order of magnitude estimate of the shear viscosity, 
it is usually correct to pick up as a trial function the case considered for ly = 0. It usually 
gives also the correct parametric behavior of the viscosity. However, in this case this is not 
so, as some values of the parameter u cause a zero-mode in the collision term^. We have 
identified two of these zero-modes. 

Let us first see why these zero modes occur. For a general trial function, the denom- 
inator of Eq. (5.19), and as explained at length in the previous sections, is dominated by 
contributions to the integral coming from (almost) collinear collisions due to the divergence 
in the naked t-channel amplitude. This divergence is however regulated by the inclusion 
of Landau damping, see Eq. (4.5). For a perfect collinear collision, the conservation of 
energy momentum at the vertex (5'^)(k -|- q — p)5{Ej^. -\- uj — Ep), where io and q are the 
energy and momentum transfer, together with Ep = vp requires that p, k and q are all 
parallel. Therefore all entering and outgoing boson momenta in a collinear collision are 
parallel. However, in solving the transport equation, a special choice of trial wavefunction 
may force a zero in this collinear limit and it can thus achieve a cancellation. Then this 
zero mode dominates the inversion of the collision operator C when solving the Boltzmann 
equation. 

This property of the (linearized) collision operator stems from the term 

Y^^iPvPn^iPvPl- (5.22) 

Now observe the cancellation for u = —2. In this case, Pijp~'^ = piPj — ^ depends only 
on the unit vector along the momentum, and 

A(PijP^) = piPj + kikj - p'ip'j - k'ik'j = . (5.23) 

for a perfect collinear collision. A second zero-mode is given by the variational value 
u = —1. Then 

MPijP") = - i EpPiPj + Ekhkj - Epip\p'j - Ek'k'ik'j + -Sij{Ep + Ek - Ep, - Ef,,) \ . 

(5.24) 
This also vanishes in the collinear limit upon employing the energy-conservation equation. 

The collision operator does not vanish because contributions to the integral not in the 
perfect collinear limit still provide a contribution, and thus a near-zero mode, not an exact 
zero mode, appears. 

We have carried out a variational study of r/[z^] numerically (see Appendix B), and we 
have reached to the conclusion that rj is maximized for the value v = —1, which corresponds 
to one of the (near) zero-modes explained above. 

We quote here the final outcome of these analysis. Due to the existence of the (near) 
zero mode with u = —1 that suppresses collinear scattering, all 2 <-^ 2 channels are included, 

^We thank the referee for calling this important point to our attention. 
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and the viscosity behaves as 



7] = 1.3- 10-^-;^ MeV^ . (5.25) 



For illustration we also give the result for low momentum transfer that is amenable to 
analytical treatment when only the t-channel is included the viscosity behaves instead as 

6. Contribution of the in-medium electromagnetism to the viscosity 

In considering transport coefficients for the CFL quark star we should examine the contri- 
bution of the in-medium electromagnetism, since the very light electron could conceivably 
transport momentum efficiently. The effective coupling of the in-medium photon and the 
electron is e = ecos^, where e is the electromagnetic coupling constant, and the mixing 
angle is [29] 

cos9= , ^^ =, (6.1) 

V3ff2 + 4e2 

CFL quark matter is a dielectric at low energies, with dielectric constant [29] 

2 f^iP 



In particular, this implies that the in-medium photon travels at a speed v-y = 1/ve^, less 
than the speed of light. Reflection and refraction properties of light on CFL quark matter 
have been studied in [30]. 

The value of the shear viscosity can be borrowed from QED, by taking into account 
the following replacements in the electromagnetic fine structure constant 

6 COS 

a ^ a = -= = a — -=— (6.3) 

47rve ve 

Transport properties in a photon-electron system in the regime T < m-e, where m-e = 
0.5 MeV is the electron mass, have been computed in [31]. While the computation, based 
on electron-in medium photon scattering at an arbitrary temperature can only be done 
numerically, approximate analytical results give 

V'-"^- = 1^- , (6.4) 

where ax is the Thomson cross-section 

and Xy^e = ^7,e/(^7 + ne) measures the relative concentrations of photons/electrons. 
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For an order of magnitude estimate, one can take a ~ 1/137, as for moderates densities 
€ never becomes too large. Then, one can see that the electromagnetic contribution to the 
shear viscosity is negligible compared to the phonon contribution at temperatures of order 
0.1 MeV. Further lowering the temperature makes Eq. (6.4) increase rapidly because of 
the Boltzmann density factor e"™'^ in the denominator containing the electron density. 
However, light by light scattering processes come to play (these can be treated with the 
Euler-Heisenberg effective lagrangian as reviewed in [32]) and the photon scattering cross 
section is not exponentially low. 

7. Further improvements 

Let us comment on possible corrections to our results. First, one could better determine 
the self-interactions of the phonons. According to ref. [19], corrections to the effective 
Lagrangian Eq. (2.8) can be obtained with a better determination of the EOS of CFL 
quark matter. Nevertheless, as the EOS is dominated by a term going as ~ ^u^, corrections 
to the phonon Lagrangian will not modify the leading parametric dependence on fj, of the 
three-phonon vertex, the relevant one for transport in the superfluid, and would only alter 
the numerical factor in front of this term in the effective Lagrangian. 

Second, one could better determine the collision term entering the transport equation. 
We saw that Landau damping regulates an otherwise divergent scattering matrix of binary 
collisions (see Sec. 3). Taking into account the width of the exchanged phonon amounts 
to a resummation of diagrams, where multiple splitting/joining scatterings are treated as 
independent classical events. Since the mean free path associated to the 2 ^^ 2 and the 1 ■s-^ 
2 collisions are of the same order, this approximation is not quite correct. One should take 
into account the Landau-Pomeranchuk-Migdal (LPM) effect, which would better determine 
the collision term [33, 34, 35]. This correction would not modify the parametric behavior 
of 7], but would again affect the preceding numerical factor.^ 

One should as well be concerned about possible changes in the shear viscosity when 
either higher derivative interactions or higher order corrections are included in the phonon 
dispersion relation. As already discussed in Sec. 2, Zarembo [20] has shown that to all 
orders in — -, the interaction terms are suppressed by powers of fi, not A. Also of interest 
to us is Zarembo's observation that the phonon dispersion relation is given by Eq. (2.9). 
As k is increased, the phonons move slower. The negative coefficient in the correction to 
the linear law above reveals that corrections to the dispersion relation beyond Son's theory 
suppress collinear splitting (a phonon cannot decay into two phonons of larger joint energy 
at T = 0). Therefore we are safe in ignoring the 1^2 processes in our computation of 
the shear viscosity. 

We note that for superfluid He^, the phonon dispersion relation used in ref. [21] to 
match the experimental value of the shear viscosity is of the form 

u;s{k) = vk[l + g{k)] . (7.1) 



We thank L. YafFe for this observation. 
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with g{k) a positive function that tends to zero for k ^ 0. In this way, phonons with high 
momenta move faster than those of low momenta and can decay into two slower phonons. 
Maris used the experimental value of the viscosity to gain information about the function 
g{k). It then turns out that in superfluid He'* almost collinear processes are more relevant 
than large-angle collisions for the shear viscosity at very low temperature. 

For astrophysical applications an order of magnitude estimate for the values of the 
viscosities is good enough to study properties of compact stars. Other transport properties 
would also be needed, but they will be the subject of a different project. 



8. Discussion 

We finally summarize our main findings and conclusions. CFL quark matter at very low T 
behaves as a superfluid of the sort He*, rather than of the sort of He^, as wrongly assumed 
in the literature. Dissipative processes in such a cold regime are dominated by phonon- 
phonon scattering. This analogy is reflected in T-power laws for the thermal properties of 
phonons of the CFL and He* superfluids. 

For astrophysical applications, we should emphasize the following. At sufficiently low 
T the phonon mean free path would exceed the radius of a compact star. We can give a 
crude estimate of the temperature when this will occur, simply by considering the equation 

i? < L ~ //T^ . (8.1) 

If the quark chemical potential is of order /i ~ 500 MeV, and we consider i? ~ 10 km, we 
find that for T < 0.06 MeV superfluid phonons do not scatter within the star. Transport 
coefficients could then be dominated by the tiny contribution of the in-medium electromag- 
netism, but an evaluation of the photon mean free path also shows that for T ~ 0.02 MeV it 
also exceeds the radius of the star [18]. Below that temperature, CFL quark matter in the 
star would behave as a perfect superfluid, as a hydrodynamical description of the phonon 
and electron fluids would be meaningless. The superfluid hydrodynamical equations would 
then be given by Eqs. (2.3) and (2.5), showing then no dissipation. 

In a rotating superfluid there are vortices. To study the rotational properties of a 
hypothetical CFL quark star, one cannot obviate that fact. In view of our results, the 
analysis of r-mode instabilities of a CFL quark star should then be redone, taking into 
account both the temperature regime of the star, and the vortex dynamics of the CFL 
phase. The bulk viscosity should also be needed, but we leave this computation for a 
different project. 
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A. Evaluation of the 2^2 processes 

In this Appendix we give some details on the Montecarlo numerical evaluations of the 
multiple integrals appearing in two-body elastic scattering. Here we treat Eq. (4.3) at 
length. The integral is dominated by t-channel scattering (the u-channel contribution is 
equal and amounts to a factor 2) due to the collinear singularity. Then for production 
runs it is sufficient to maintain the t-channel scattering matrix |Tp and neglect the rest of 
the interaction diagrams. The numerator of Eq. (4.3) has twelve integrations. Three, for 
example over k', can be immediately performed with the help of the three-momentum 6 in 
Eq. (4.3). We choose for the nine remaining integrals the incoming three- momenta p, k, 
and the transferred momentum q = p — p' = k' — k. 

Our freedom to choose the third axis makes the angular integrals around q to be An. 
Since we approximate the integral to be dominated by t-channel exchange, the scattering 
matrix element |Tp can only depend on the two polar angles (through the derivative 
coupling in the vertex), p-q and k-q and not on the azimuthal angles (pp and (p^^. Therefore 
the latter are also trivial. 

Introduction of an auxiliary variable to (energy transfer) allows to write the energy 5 
in Eq. (4.3) as 

/ + 00 
doj S{w -Ep + Ep,)5{uj + Ek-Ek') . (A.l) 

-oo 

Then it is easy to employ the resulting two 6 functions to compute the remaining polar 
angular integrals. In exchange, the integral over the auxiliary variable remains. With these 
manipulations we obtain 

1 />oo /-vq /"OO /•— oo 

r(2^) = ^3 3 r^W8 dq du; dk dpUUl + Mil + /pOirp , (A.2) 

where |Tp is given in Eq. (4.5). This is an integral in four dimensions dominated by 
the region around the on-shell singularities for lo c^ vq and lj ~ —vq. To obtain good 
numerical convergence these singularities, which are regulated by Landau damping, need 
to be isolated. First, perform one more change of variables to 

X = uj + vq ; y = to — vq . (A. 3) 

The integrals over x and y run from to oo and from — oo to respectively. The on-shell 
singularities occur now at x = 0, y = 0. Therefore we introduce a low momentum cut-off 
on these variables. A, that serves to split the regions of integration in the xy plane as (1) 
(0,A) X (-A,0); (2) (0,A) x (-oo,-A); (3) (A, oo) x (-A,0); (4) (A, oo) x (-oo,-A). The 
value of the integral should of course be independent of the value of A since this merely 
splits the integration region. Choosing A ^ T allows to evaluate most of the integrand in 
regions (2) and (3) at x = , y = respectively, and therefore the x (y) integral can be 
performed analytically. Explicit numerical evaluation shows that regions (1) and (4) are 
negligible with respect to the contributions from regions (2) and (3) whenever A is large 
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imni 



compared to — ^^^- (or the same expression reversing x and y). Under these conditions 
we need to concentrate only on regions (2) and (3). 
For region (3) the relevant y integral is 

.0 1 
I^F{P,x,y = 0) dy ^ —- . (A.4) 

Since in this region x > A, the imaginary part takes a constant value as y ^ 0. Therefore 
the integral is regulated and yields an arctangent. All the derivative couplings and Bose- 
Einstein factors are smooth at y = and have been pulled out of the integral as a constant 
(still function of x) F(P,x,y = 0), since the interval is small respect to T. Therefore 

/ = F{P,x,y = 0) ^ j/^ . (A.5) 

Finally the damping rate can be given to very good accuracy in terms of three quadratures 

^ (fs+fg) (A.6) 



128?(3)t;2T3 (Slfi^y 

/"OO /"OO f'OO l<r\ 

T3^ dx dp dkfphil + /kO(l + /pO , „/ -r {F{P, -Q)F{K, Q)) |,=o • 

Jx Jx/2v Jo xlmn(a;,y = 0) 

The remaining three integrals are routinely calculated with a computer code and yield the 
damping rate. Evaluation of the shear viscosity is dealt with in the next two sections. 

B. Numerical evaluation of shear viscosity 

Our Montecarlo program, employing the subroutine Vegas [36] has been written in two 
versions. One perfoms three and four dimensional integrals keeping only the t-channel 
exchange contribution to the cross section as previously described, and is useful to isolate 
this precisely. If other scattering processes also contribute, this program cannot be used 
since two azimutal angles are now not trivial. In this case a second version performs the 
five dimensional integral over all phase space and is a generic purpose integrator (but 
diluting the available computing power over a larger region makes it less precise to capture 
t-channel singularities). In table 2 where we present various numerical results, these two 
versions are denoted as (t-channel) and (all) respectively. 

The numerical results from our Montecarlo for the shear viscosity are given in full in 
table 2. The variational calculation is somewhat complex and care needs to be exercised 
upon its interpretation. In this table, the first column gives the temperature in MeV. The 
second is the variational parameter z^ that varies from a strong infrared enhancement to 
polynomial suppression. The fourth column displays the values obtained for the viscosity 
with all interactions included. These are extracted in the third column. 

The fourth column displays very clearly the ti^ behaviour characteristic of the naive 
power counting. However, for a general i', one should go to the result of the t-channel 
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Table 1: Numerical cancellation of the internal A dependence in ?](— 2) given by eq. (5.19) em- 
ploying the zero mode trial wavefunction v = —2. The temperature is T = 0.001 McV, the cutoff 
for external momenta is A = lOOT and the cutoff over the x, y integrations is 0.6A, large enough 
to cover all relevant phase space but small enough to avoid the Montecarlo points to be dispersed 
where the Bose-Einstein factor damps the integrand. The calculation is performed with 256000 
points in each integration domain, doubled four times for convergence, and the integral is com- 
puted ten times for each fixed number of points. The standard deviation obtained is 1% . A small 
logarithmic drift of the result is still observed, but this is now of the order of the uncertainty in 
the Montecarlo and not very relevant. A meaningful cutoff A needs to be smaller than T but of a 
similar order of magnitude. 



A 


7]{U 


= -2){MeV^) 


T/12 




1.09 


1033 


T/25 




1.10 


1033 


T/50 




1.12 


1033 


T/lOO 




1.14 


1033 


r/200 




1.16 


1033 



exchange. For example, for v = Q that corresponds to the usual zero'th order approximation 

4 

in a polynomial expansion, and at T = 0.1 MeV, the viscosity would behave as ry ~ 0.8-^. 

But this conclusion would miss the effects of the zero mode. By studying r]{u) with 
the t-channel calculation (again the third column, at fixed T), one observes that the zero 
modes discussed above corresponding to z^ = —2 and v = —1 yield variational maxima, 
especially u = —1, a, mildly dependent on the numerical grid and cutoffs. The analytical 
considerations suggest a logarithmic dependence, as discussed in the main text. 

Finally a glance at the fourth column reveals that this formula is much larger (and 
therefore the integral (5.21) very suppressed) than the corresponding naive counting for 
this value of the variational parameter. This entry also happens to be the maximum of 
the third column, finally concluding that the viscosity must behave in this temperature 
range as in Eq.(5.25). The plot 4 enlightens this discussion by showing how the t-channel 
contribution is much lower than the corresponding T~^ behaviour for v = —1, u = —2. 
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Figure 4: The variational evaluation including boson exchanges in the interaction shows a max- 
imum corresponding to the zero modes that make the t-channcl contribution suppressed. The 
maximum within the functional family - — ^^ ^^ corresponds to ;/ = — 1. The Temperature is fixed 
at 0.01 MeV. 
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